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Introduction: The initial accretion of primitive bodies 
(here, asteroids in particular) from freely-floating nebula par- 
ticles remains problematic. Traditional growth-by-sticking 
models encounter a formidable “meter-size barrier" [1] (or 
even a mm-to-cm-size barrier [2]) in turbulent nebulae, making 
the preconditions for so-called “streaming instabilities" diffi- 
cult to achieve even for so-called “lucky" particles [3]. Even 
if growth by sticking could somehow breach the meter size 
barrier, turbulent nebulae present further obstacles through the 
l-10km size range [4], On the other hand, nonturbulent nebu- 
lae form large asteroids too quickly to explain long spreads in 
formation times, or the dearth of melted asteroids [5]. Theoret- 
ical understanding of nebula turbulence is itself in flux; recent 
models of MRI (magnetically-driven) turbulence favor low- 
or-no-turbulence environments [6], but purely hydrodynamic 
turbulence is making a comeback, with two recently discov- 
ered mechanisms generating robust turbulence which do not 
rely on magnetic fields at all [7,8]. 

An important clue regarding planetesimal formation is 
an apparent 100km diameter peak in the pre-depletion, pre- 
erosion mass distribution of asteroids [9]; scenarios leading 
directly from independent nebula particulates to large objects 
of this size, which avoid the problematic m-km size range, 
could be called “leapfrog" scenarios [10-12]. The leapfrog 
scenario we have studied in detail involves formation of dense 
clumps of aerodynamically selected, typically mm-size parti- 
cles in turbulence, which can under certain conditions shrink 
inexorably on 100-1000 orbit timescales and form 10- 100km 
diameter sandpile planetesimals. There is evidence that at least 
the ordinary chondrite parent bodies were initially composed 
entirely of a homogeneous mix of such particles [4], Thus, 
while they are arcane, turbulent concentration models acting 
directly on chondrule size particles are worthy of deeper study. 
The typical sizes of planetesimals and the rate of their forma- 
tion [11,12] can be estimated using a statistical model with 
properties inferred from large numerical simulations of turbu- 
lence [13], Nebula turbulence is described by its Reynolds 
number Re = ( L/r) ) 4//3 , where L = Ha 1 ^ 2 is the largest 
eddy scale, H is the nebula gas vertical scale height, a the 
turbulent viscosity parameter, and p is the Kolmogorov or 
smallest scale in turbulence (typically about 1km), with eddy 
turnover time t v . In the nebula, Re is far larger than any 
numerical simulation can handle, so some physical arguments 
are needed to extend the results of numerical simulations to 
nebula conditions. In this paper, we report new physics to be 
incorporated into our statistical models. 

Cascade model and multiplier distributions: The spa- 
tial distribution of particle concentration can be captured sta- 
tistically by a cascade model [11,13] which predicts the proba- 
bility distribution functions (PDFs) for dense particle clumps; 
these PDFs are essentially the volume fractions of the nebula 
which have the necessary properties (solids mass and local 
vorticity) for planetesimal formation. A cascade model pre- 
sumes that, as energy flows from large eddies to smaller ones, 
particles and fluid properties are partitioned unequally at each 
“level" of the cascade from “parent" into “daughter" eddies. 
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Figure 1: The symbols indicate the multiplier pdfs 

for turbulent dissipation in high- .Re atmospheric turbu- 
lence [10], showing scale invariance over 300-3000?7. The 
red /3-functions (/3=3, 6, and 48) provide context for /3 
values in figure 2. In [11] we adopted level-independent 
multiplier pdfs with /3 ~ 3 based on our own simula- 
tions. 

The probability distribution functions (pdfs) of the so-called 
multipliers by which particle and fluid properties are parti- 
tioned have widths that are parameterized by the parameter 
/3; smaller /3 means larger width ([13]; see figure 1 above), 
or a higher probability of strongly asymmetrical partitioning 
which, repeated over a number of levels, leads to clumpier or 
more intermittent local values of particle and fluid properties. 

The wide range of scales between L and p is loosely 
called the inertial range; in this range the equations of motion 
are scale-free, and the multiplier pdfs for dissipation of tur- 
bulent kinetic energy (which, like the chondrule-size particles 
modeled by [11-13], is characterized by Kolmogorov eddy 
timescales) are known to be scale-independent in atmospheric 
flows from hundreds to thousands of Kolmogorov lengthscales 
[14;figure 1], This fact motivated [1 1,13] to develop acascade 
model in which the multiplier pdfs for particle concentration 
<f> were also independent of turbulent eddy lengthscale. Subse- 
quent disagreement with results of others [15] led us to a more 
in-depth study using highly resolved fluid simulations [16,17] 
that revealed a scale dependence of multiplier pdfs, which will 
affect the formation rate and masses of primary planetesimals 
in the scenario of [11,12], The approach is described in [18]. 
Here we present further analysis along these lines. 

Results: Figure 2 shows our results for the widths of the 
multiplier pdfs [18], binned on a wide range of lengthscales r 
given in units of p. In this figure, the particle stopping time 
t a has been normalized to a scale-dependent value St r using 
the inertial range timescale ratio ( r/p)~ 2 ' 3 , and the constant 
parameter St = t s /t v where t v is the Kolmogorov eddy time. 
In this plot, a value of St r = 1 represents a particle having 
stopping time equal to the eddy time for the spatial scale r. 
Except at the largest binning scales (r > 128t/), this scaling 
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Figure 2: (3 (like inverse width) of multiplier pdfs for 
particle concentration, plotted as a function of St r 
(stopping time t s normalized by the turnover time of 
an eddy with scale r), for a variety of binning scales r 
ranging from 1277 to 512?; ~ L/ 2. With this scaling, 
the multiplier pdf (3 values collapse to a universal curve 
except at the largest scales r > L/10. The scale-free 
(3 ~ 3 for dissipation (figure 1) is shown. A universal 
curve allows our cascade models to be extended to scales 
smaller than resolvable by current numerical models. 

collapses our previous results [18] into an invariant family of 
curves, allowing us to predict the value of (3 for any combina- 
tion of lengthscale and particle stopping time. The minimum 
near St r = 0.2 suggests that particles with stopping time 
somewhat shorter than the eddy time on lengthscale r are most 
strongly concentrated within spatial averages over lengthscale 
r. This scaling collapses the multipliers much better than an 
alternate scaling found by [19] which collapses the normal- 
ized concentration pdfs themselves; the concentration pdfs 
can be thought of as the cumulative result of multipliers acting 
sequentially over all larger scales. After a number of eddy 
bifurcations to smaller scale, particles of a small but finite 
range of sizes might thus end up concentrated within a sin- 
gle clump, which may be consistent with new determinations 
of chondrule size distributions which find them to be broader 
than previously believed [20], Our model will ultimately make 
testable predictions along these lines. 

However, the curves for the largest lengthscales (r > 
L/10 where L ~ 1000 rj) diverge from this scale-invariant 
curve, probably because, in the top decade of lengthscales, 
turbulent stretching and vortex tube formation have not yet 
reached their fully developed state (see [8] and references 
therein). For these largest scales, multiplier statistics are better 
characterized by their scale as a fraction of L rather than as a 
multiple ofr/. Moreover, at the largest scales, where nominal 
multiplier (3 values are large, the /3-function itself does not pro- 
vide a complete description of the cascade properties (figure 
3). The scale at which the /3-function becomes an adequate fit 
to the actual multiplier pdf is close to the scale at which the 
f3(St r ) curves of figure 2 become universal. 

The cascade model for planetesimal formation by turbu- 
lent concentration [11] can now be modified to allow for the 
observed level-dependence in the particle concentration multi- 
plier pdfs (figures 2 and 3), as well as in the multiplier pdfs for 
gas enstrophy (vorticity squared) which also determine the cas- 
cade results and, we have found, also differ from values used 



Figure 3: Multiplier pdfs for nominally St — 1 parti- 
cles, at four different scales. The red curves are the 
best-fitting overall /3-functions, which get narrower for 
a given particle size to large scales (in figure 2, points 
move up and to the left on the universal collapse curve) . 
At smaller scales, lying within the universal collapse do- 
main, the /3-function itself provides a good (and scale- 
invariant) fit to the actual multiplier pdf. However the 
actual pdfs (black curves) have broader exponential tails 
at large spatial scales. 

in [11,12]. However, the asymptotic (3 value we determine 
for dissipation does not quite match the scale-free value deter- 
mined by [14] in atmospheric flows. Thus, some effort must 
be dedicated to understanding how the results presented here 
may change with Re itself, and we are pursuing this. Because 
nebula properties such as turbulent intensity and gas density 
are unknown to an order of magnitude [11,12], these findings 
do not, on the face of it, qualitatively change the outlook for the 
scenario, but will lead to quantitative changes in predictions of 
planetesimal IMF and chondrule size distributions. 
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